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ABSTRACT 


V 


Synthetic  data,  simulating  the  effects  of  atmospheric 
extinction  on  GEODSS  observations,  are  presented.  Approximate 
representations  for  the  atmospheric  extinction  coefficients  and 
for  spectral  photon  distributions  are  also  developed.  Together 
these  data  are  intended  for  evaluation  of  the  effects  of  the 

«-<*->>-  . . . i . t.  . W - ■ ' *4 

atmosphere  on  SOI  data,  and  for  the  development  and  testing  of 


models  for  correcting  atmospheric  extinction. 
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I. 


INTRODUCTION 


The  correction  of  astronomical  brightness  measurements  for 
the  effects  of  atmospheric  extinction  is  a complicated  process 
for  which  no  universally  satisfactory  procedure  exists.  The 
problem  is  especially  difficult  in  the  case  of  the  establishment 
of  a catalog  of  reference  standards.  Here,  mediocre  extinction 
corrections  can  lead  tc  unacceptably  large  systematic  errors  which 
destroy  the  internal  consistency  of  the  catalog.  The  problem  has 
been  extensively  studied  for  the  conditions  appropriate  to  typical 
astronomical  measurements.  The  conditions  of  GEODSS  measurements, 
however,  are  rather  different  and  the  appropriate  methodology  for 
these  corrections  is  not  well  developed.  See  references  (1)  and 
(2)  for  discussions  of  the  fundamental  correction  process. 

There  are  two  approaches  to  making  extinction  corrections: 
direct  integration  and  modeling.  If  the  spectral  distribution, 
1(A),  of  the  source  being  measured  and  the  monochromatic  extinc- 
tion coefficient,  k(A),  are  known  then  a broadband  extinction  may 
be  calculated  by  direct  integration.  If  1(A)  and  k(A)  are  among 
the  unknowns  sought,  an  iterative  technique  can  be  expected  to 
converge  given  enough  data.  Unfortunately,  the  amount  of  data 
required  for  this  to  work  rules  it  out  in  the  case  of  GEODSS  and 
even  most  astronomical  observations.  Modeling  works  by  finding 
useful  (i.e.,  some  tradeoff  between  accurate  and  simple)  approxi- 
mations to  such  quantities  as  k ( A ) and  developing  programs 


specifically  designed  to  determine  the  parameters  needed  to 
make  extinction  corrections.  The  GEODSS  program  needs  two  rather 
different  extinction  models.  One  must  be  highly  precise  and  need 
not  be  particularly  simple,  the  second  must  be  capable  of  being 
evaluated  using  very  little  observing  time  and  can  tolerate 
relatively  poorer  precision. 

This  paper  has  been  prepared  in  order  to  facilitate 
development  of  extinction  correction  models  for  GEODSS  photometry. 
Section  II  contains  values  of  broadband  extinction  calculated  by 
direct  integration.  These  are  the  values  that  any  proposed 
extinction  model  must  be  capable  of  producing.  The  model  must 
also  work  in  the  other  direction  and  be  able  to  provide  the 
relevant  parameters  on  the  basis  of  observational  data  contamina- 
ted with  noise.  To  test  this  ability  a set  of  synthetic  data  has 
been  included.  Section  III  presents  approximate  functional 
representations  of  the  relevant  physical  quantities.  These 
provide  the  form  and  initial  estimates  needed  in  the  development 
of  our  extinction  model.  Finally,  the  data  presented  here  should 
prove  useful  in  the  evaluation  of  error  budgets  for  an  operating 
system. 


II.  SYNTHETIC  OBSERVATIONS 

The  measured  magnitude  of  an  object  is  given  by 
m'.  = m . + 6 ( A ; j ) + e 

J oj 

where  m^  = measured  magnitude  through  filter  j 

mQj=  exo-atmospheric  (catalog) , j-filter  magnitude 
6 (A; j ) = atmospheric  extinction  for  the  j-filter  magnitude 
as  a function  of  airmass,  A 

e = error  - assumed  to  be  drawn  from  normal  distribution 
In  order  to  create  a set  of  m'^  it  is  first  necessary  to  calculate 
values  of  6(A;j)  and  to  generate  a set  of  random  errors.  The  6 
may  be  calculated  from 


6 (A; j ) = --In 


/ 1 (X ) T (A, X ) Rj  (X ) dX 
/ 1 (d)  Rj  ( X ) dX 


where  I ( X ) = the  spectral  distribution  of  the  object 

T (A, X ) = the  atmospheric  transmission  as  a function  of 
wavelength  and  airmass;  T(A,X)  = e qAlc^ 

A = the  airmass,  a function  of  zenith  distance 
k(X)  = the  extinction  coefficient,  in  magnitudes,  as  a 
function  of  wavelength 

Rj (X)  = the  instrumental  response  profile  for  filter  j 
q = . 4 In  10 

The  integrals  may  be  evaluated  numerically  for  known  I(X). 
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The  values  needed  to  carry  out  the  integrations  are 


tabulated  as  Table  I.  R^,  and  R^ , are  the  unfiltered  GaAs 
photocathode  response  and  the  GaA_  response  modified  by  the 
addition  of  a Schott  GG435  filter.  This  modification  was  made 
in  an  attempt  to  make  the  extinction  characteristics  more 
tractable.  R^  and  R^  are  the  effective  profiles  of  the  W-  and 
X-filters  of  the  Lincoln  three-color  system.  These  are  defined 
by  1mm  BG18  plus  1mm  BG27  for  W and  3mm  OG590  plus  2mm  KG3  for 
X,  with  a GaAs  photocathode.  The  central  wavelength  and  full 
width  at  half  power  are  respectively  43008  and  14008,  and 
6720&  and  120o8.  Five  spectral  distributions  are  listed: 
sunlight,  sunlight  reflected  from  gold,  and  black  body  radiation 
at  4400°K,  5800°K,  and  6900°K.  These  temperatures  correspond  to 
the  effective  temperatures  of  stars  of  spectral  types  K2,  G2 , and 

: 

F2,  i.e.,  the  sun's  type  and  one  type  to  either  side.  The 
atmospheric  extinction  coefficients  are  for  a "typical" 
atmosphere  as  described  in  appendix  B of  reference  (3) . 

Values  of  the  ratio 

f I (A  JR..  (X ) T (A , A ) dX 
/ I ( X ) R j (d) dX 

were  calculated  for  several  values  of  the  air  mass  and  for  all 
combinations  of  1(A)  and  R j ( X ) . These  values  were  converted  into 
astronomical  magnitudes  and  are  tabulated  as  Table  II.  The  first 
section  of  the  table  lists  values  of  the  atmospheric  extinction 
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TABLE  II 


VALUES  OF  THE  ATMOSPHERIC  EXTINCTION  AS  A FUNCTION  OF 
AIR  MASS  FOR  A VARIETY  OF  SOURCE  SPECTRAL  DISTRIBUTIONS. 

For  sunlight,  the  column  labeled  X'  contains  values  calculated 
with  the  "extra"  extinction  at  X = . 76y  removed. 


a.  SUNLIGHT 


Air  Mass 

E 

E' 

W 

X 

X' 

1.  0 

0™372 

0™288 

0?555 

0m252 

0I?238 

1.25 

.456 

. 358 

.690 

.315 

.297 

1.5 

. 538 

.428 

.824 

. 377 

. 356 

1.75 

.617 

.497 

.956 

.439 

.415 

2.0 

.695 

. 565 

1.088 

. 500 

.474 

2.5 

.847 

.701 

1.347 

.622 

.591 

3.0 

.993 

. 834 

1.601 

.744 

.708 

3.5 

1.135 

.965 

1.852 

.864 

.825 

b.  : 

SUNLIGHT  - 

GOLD 

E 

E' 

W 

X 

1.0 

0™309 

0^268 

0™524 

0^251 

1.25 

.381 

. 334 

.652 

. 313 

1.5 

.451 

.399 

.778 

.375 

1.  75 

.521 

.464 

.903 

.436 

2.0 

.589 

.528 

1.027 

.498 

2.5 

.724 

.655 

1.271 

.619 

3.0 

.855 

.780 

1.511 

.740 

3.5 

.983 

.903 

1.748 

.860 
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TABLE  II  (Continued) 


c.  BLACK 

- BODY 

4400°K 

Air  Mass 

E 

E' 

W 

X 

1.0 

0™315 

0™266 

0?541 

0?250 

1.  25 

.388 

.331 

.673 

. 312 

1.  5 

.460 

. 395 

.804 

.373 

1.75 

. 530 

.459 

.933 

.434 

2.0 

.599 

. 523 

1.061 

.495 

2.5 

. 734 

.648 

1.313 

.615 

3.0 

.865 

.772 

1.561 

.735 

3.  5 

.993 

.893 

1.805 

.854 

d . BLACK 

- BODY 

5800°K 

E 

E' 

W 

X 

1.0 

0?395 

0™286 

0?574 

0?252 

1.25 

.482 

.356 

.714 

.314 

1.  5 

.566 

.425 

.851 

. 376 

1.  75 

. 647 

.494 

.988 

.438 

2.0 

.727 

.562 

1.122 

.499 

2.5 

.881 

.697 

1.388 

.621 

3.0 

1.029 

.829 

1.649 

.742 

3.5 

1.172 

.960 

1.905 

.863 
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TABLE  II  (Continued) 


e . BLACK 

- BODY 

6900°K 

Air  Mass 

E 

E' 

W 

X 

1.0 

0I?452 

0?295 

ol^i 

0^25  3 

1.25 

.547 

. 368 

.736 

. 316 

1.5 

.639 

.439 

.878 

.378 

1.75 

.727 

.510 

1.018 

.440 

2.0 

.813 

.580 

1.157 

.502 

2.5 

.979 

.719 

1.429 

.625 

3.0 

1.137 

.856 

1.696 

.746 

3.5 

1.289 

.990 

1.958 

.867 
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for  sunlight  as  seen  by  the  four  response  profiles  R^, , RE,,  R^, 
and  R^,  and  an  additional  column  labled  X' . This  column  contains 
values  of  extinction  arising  from  a modified  atmospheric  extinc- 
tion - one  having  the  "extra"  extinction  due  to  the  atmospheric 
A-band  removed.  This  has  been  done  to  facilitate  later  compari- 
son with  a simple  model  which  does  not  include  the  effect  of 
the  A-band.  Two  effects  of  the  wavelength  dependence  of 
atmospheric  extinction  are  immediately  evident.  First,  the 
extinction  in  the  bluest  filter  (W)  is  much  the  highest,  and 
second,  for  any  response  bluer  objects  (here,  black  body 
radiation  at  higher  T)  suffer  greater  extinction. 

Random  errors  were  generated  in  the  following  way.  A 
pseudo-random  number  generator  was  programmed  to  provide  numbers 
evenly  distributed  over  the  interval  -.5  < y < + .5.  A value  of 
t was  found  such  that  |y|  = (2-n)~^  /^e  x ^2dx,  and  the  sign  of  y 
was  then  assigned  to  t.  The  numbers  t are  thus  drawn  from  a 
normal  distribution  with  a mean  of  zero  and  standard  deviation 
of  unity.  In  fact,  for  the  values  of  t used,  t = -0.007  and 
a = 1.059.  Values  of  t were  then  multiplied  by  a typical  mean 
error  of  a single  observation  and  added  to  the  magnitude  to 
produce  the  synthetic  data.  Two  typical  values  were  used,  a 
constant  0™030  and  an  air  mass  proportional  0™010  A^ . 

Synthetic  data  were  generated  for  ten  "objects":  two  black- 
body  sources  at  each  temperature,  three  sunlight  sources,  and  a 
sunlight  from  gold  source.  Instead  of  using  the  two  magnitudes 
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TABLE  III 


CATALOG  MAGNITUDES  FOR  OBJECTS  FOR 
WHICH  SYNTHETIC  OBSERVATIONS  WERE  GENERATED 


OBJ 

m 

e,x  o 

c 

o 

REMARKS 

1 

6.300 

0™715 

BB,  T = 4400° K 

2 

7.250 

0.715 

BB,  T = 4400°K 

3 

4.650 

0.064 

BB,  T = 5800°K 

4 

5.900 

0.064 

BB,  T = 5800°K 

5 

5.550 

-0.269 

BB,  T = 6900°K 

6 

6.850 

-0.269 

BB,  T = 6900°K 

7 

5.800 

0.101 

Sunlight 

8 

6.350 

0.101 

Sunlight 

9 

7.150 

0.101 

Sunlight 

10 

6.000 

0.944 

Sunlight  reflected 
from  gold 

and  m^,  the  usual  practice  of  using  a magnitude,  m^,  and  a 
color,  C = m^  - m^,  has  been  followed.  For  simplicity  the  X- 
magnitude  and  the  E-magnitudes  are  taken  to  be  equal.  The 
catalog  values  for  these  ten  objects  are  listed  in  Table  III.  The 
synthesized  observations  are  listed  in  Table  IV.  Data  are  organ- 
ized into  three  "runs"  for  each  object.  Runs  1 and  2 use  an 
irregular  selection  of  five  values  of  airmass  spanning  the  range 
of  values  calculated.  Runs  3 all  use  the  same  regular  series 
of  six  airmasses.  For  runs  1,  the  constant  standard  deviation 
error  was  added,  while  for  runs  2 and  3 the  airmass  proportional 
error  was  applied. 


SYNTHETIC  OBSERVATIONS  OF  TEN  OBJECTS 


TABLE  IV  (Continued) 


TABLE  IV  (Continued) 


TABLE  IV  (Continued) 


Ill . 


APPROXIMATE  REPRESENTATIONS 


It  is  useful  to  have  functional  representations  for  some  of 
the  data  in  Table  I.  Presented  here  are  simple  functions  which 
approximate  the  atmospheric  extinction  coefficient,  the  solar 
photon  flux,  and  the  black-body  photon  flux  for  T = 5800°K.  These 
representations  can  be  used  for  such  purposes  as  developing  models 
of  extinction  and  evaluating  specific  effects  of  atmospheric  ex- 
tinction on  real  data. 

For  the  case  of  the  extinction  coefficient,  simple  polynomi- 
als of  order  three  or  less  do  a rather  poor  job  of  approximating 
the  wavelength  dependence.  A good  fit,  however,  is  obtained  from 

k (A  ) = 0 . 096/  (A  - .26) 

Figure  1 shows  the  fit.  The  A-band  (A  = .76y)  is,  of  course, 
completely  missed  by  this  representation.  There  is  precedent  for 
this  form;  astronomers  have  long  used  a simple  1/A  extinction 
law  in  studies  of  the  interstellar  medium.  Generalizing  this  gives 

k ( A ) a — — + 6 

A - A’ 

Varying  a, A ' , and  8 should  provide  a good  approximation  to  a wide 
range  of  atmospheric  conditions.  Treating  A'  as  a free  parameter, 
however,  introduces  enormous  complications  into  the  analysis.  It 
also  seems  likely  that  the  origin  of  the  value  of  A ' lies  in 
microscopic  rather  than  macroscopic  processes  and  would  thus  not 
be  expected  to  vary  greatly.  Therefore,  it  seems  reasonable  to 
set  A'  = . 26y  and  also  to  expect  a > 8 > 0. 
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approximation  to  the 


Using  this  representation  for  k(A),  it  is  possible  to  write 


a 


s* 


power  series  expansion  for  the  atmospheric  transmission, 
transmission  law  is 

T(A,X)  = e-qAK(>) 

when  k(A)  is  expressed  in  astronomical  magnitudes  per  air 
T (A, A)  can  be  written  in  the  form 


T (A,  A ) = T(A,Aq)  [l  + ax(A)  (— ^ Ao)  + a (A)  (— ~-io)  2 + ... 

o 2 o 

Evaluating  the  coefficients  gives 


T(A,Aq)  = e_qA(Ko 

al  = AoqAKo 
a2  = V\*[(qA)2K2 

a3  * j:AoC'qA)*K= 
a4  ' J:A;UqA)'K; 


+ B) 

2(qA)Ko] 

- 6 (qA) 2K2  + 

- 12  (qA)  3K3q  + 


6(qA)Ko] 

36 ( q A ) 2 K 2 
o 


24  (qA)Kj 


where 


A 


o 


K 

o 


A - ,26y 

o 

a 

A - . 26y 
o 
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The 


mass . 


The  solar  photon  flux  may  be  fairly  well  approximated  by 
either  a cubic  polynomial  or  a black-body  distribution.  Figure  2 
shows  measured  values  of  the  solar  flux  compared  with  a black-body 
distribution  at  T = 5800°K  and  the  cubic  polynomial, 

1(A)  = -53.23(1.00  - 5.12A  + 7.11A2  - 3.18A3), 

calculated  for  the  range  .34  < A < .86.  Note  how  badly  it  fails 
outside  this  range.  It  is  useful  to  express  the  polynomial  in 


t ' 

* 

powers  of  (A  - X q ) / X q 

instead  of 

X . For  X = 
o 

. 55y 

, for  example, 

1 « 

I = 10 . 34  [l . 00  + . 

52  (A  yA-^)  - 

2 . 90  ( A -~v^)2 

+ 2. 

72  (A  ~xAQ)  3] 

o 

o 

o 

As  an  additional 

comparison. 

a cubic  was 

fit 

to  T = 5800°K 

black-body  radiation  in  the  range  .34  < A < .86.  The  function 
obtained  is 

I = 3.145(1.00  - 5.31A  + 6.68A2  - 2.64A3) 
or,  again,  setting  Aq  = . 55y 

I = 1.07(1.00  + . 58  (X  ~ X°)  - 2. 07(.-~--°)2  + 1.30(AlA°)3] 

A o A o A o 

Figure  3 shows  this  polynomial  fit  compared  to  the  black-body  curve. 

Again,  the  cubic  deviates  from  the  correct  values  outside  the 
valid  range,  although  not  as  badly  as  in  the  case  of  the  fit  to 
the  other  solar  flux  in  Figure  2.  Not  surprisingly,  the 
coefficients  for  the  black-body  fit  are  similar  to  those  for  the 
solar  fit  except  for  the  (inconsequential)  scaling  factor. 

I I 
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Figure  3.  Black-body  Radiation  at  5800 


w • - 

' 

Finally,  the  actual  black-body  functions  can  be  expanded  as 
an  infinite  order  polynomial 

B(A,T)  ~ 1 + b1  (T)-^-~_A.o+  b2  (T)  (~x~-X-°) 2 + ... 

o o 

Values  of  the  first  four  coefficients  are  given  in  Table  VI  for 
T = 6900°K,  5800°K,  and  4400°K,  and  for  Xq  = .43u,  .55u,  and  .67y. 
The  central  compartment  of  Table  V,  containing  the  coefficients 

, for  X = . 55y  and  T = 5800°K,  should  be  compared  to  the  coeffi- 

r 0 

cients  of  the  cubic  polynomial  representations. 

At  this  point  the  meaning  and  intended  use  of  the  expansions 
for  the  flux  must  be  examined  more  carefully.  The  expansions  will 
generally  be  used  to  evaluate  integrals  in  the  form 
/I ( A ) f (A ) d A = Fb  F , 

where 

F = / (A  A.P.)nf  (A)dA 
n Ao 

The  bn  are,  of  course,  just  the  coefficients  in  the  power  series 
expansion  of  I, 


o 


and  the  can  be  expected  to  become  very  small  as  n increases. 
Consider  now  the  case  of  5800°K  black-body  radiation  and  its  two 
expansions:  the  cubic  fit  defined  by  the  three  coefficients  de- 
rived above,  and  the  power  series  defined  by  an  infinite  number  of 


1 
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TABLE  V 


COEFFICIENTS  FOR  THE  POWER  SERIES  EXPANSIONS 
OF  THE  BLACK- BODY  PHOTON  DISTRIBUTION. 

X 


1 T 

. 43y 

. 55y 

.67u 

4400°K 

» 

\ 

3.61 

1.96 

.92 

b2 

.92 

-1.99 

-2.41 

b3 

-6.08 

-1.80 

1.08 

b4 

2.35 

4.87 

2.19 

' 5800°K 

bl 

1.79 

.56 

- .20 

b2 

-2.14 

-2.29 

-1.60 

•1 

b3 

-1.28 

1.83 

2.67 

b4 

4.61 

0.86 

-1.92 

6900°K 

bl 

00 

v£> 

- .12 

- .74 

b2 

-2.40 

-1.70 

- .75 

b3 

1.15 

2.64 

2.41 

2.08 

-1.65 

-3.13 

coefficients,  the  first  four  of  which  are  shown  in  Table  V • 


Now 


within  the  range  .34  i H .86  the  three  cubic  coefficients  give 
a better  fit  to  the  data  than  the  first  three  coefficients  from 
Table  V • Nevertheless,  a better  representation  of  integrals  such 
as  the  one  shown  above  would  be  obtained  by  using  the  values  from 
Table  V . Thus,  even  if  only  the  first  few  coefficients  make  a 
significant  contribution  to  the  series  representation  of  the  in- 
tegral, they  should  be  evaluated  via  a higher  order  (ideally,  in- 
finite order)  fit  to  the  data. 

When  1(A)  is  unknown,  two  additional  problems  arise.  First 
1(A)  will  be  characterized  or  characterizable  by  an  intensity  and 
a few  color  indices.  In  order  to  evaluate  the  integrals,  it  will 
be  necessary  to  relate  the  bn  to  these  color  indices.  This  leads 
naturally  to  the  next  problem  - the  non-independence  of  the  bn . 

The  analysis  used  to  determine  1(A)  is  based  on  the  assumption  of 
independent  parameters.  When  the  parameters  are  not  independent 
the  values  derived  for  them  will  be  unreliable.  In  the  case  of 
black-body  radiation,  there  is  only  one  free  parameter,  T (or 
equivalently  b^.  For  real  radiation,  the  strict  functional  de- 
pendence on  only  one  parameter  no  longer  holds,  but  neither  is 
there  a set  of  many  truly  independent  bn.  It  is  generally  pre- 
ferable to  accept  a slightly  poorer  fit  to  the  data  using  a minimum 
of  parameters  than  to  improve  the  fit  by  adding  free  parameters  of 
questionable  nature. 
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